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consistent across different test methods and are also in good agreement with the results obtained by 
others. A summary of engineering properties of Sydney Hawkesbury sandstone is given in Table 


4, 
Tabdie +. Engineering properties of Sydney Hawkesbury sandstone 
Parameters | Testmettods | pico pa Range of values Average value 
| | asad | 8.7 ~23.1 (MPa) 16.9 (MPa) 
o, | uniaxial airdned | 25.1 ~356 (Pa) | 30.1 (MPa 
[ _ovendried _| 8.1 ~ 46.2 (MPa) 41.3 (MPa) 
high humidity ine - 
ea 2.18 ~ 6.45 (GPa) 4.50 (GPa) 
E uniaxial | airdried [5.30 ~ 8.60 (GPa) 6.79 (GPa) 
| _ovendried | 7.40 ~ 10.2 (GPa) 8.42 (GPa 
{ov triaxial “saturated | 0.09~0.10 | 0.10 
Srinsacey triaxial | saturated | 48.5 ~ 50.0 (deg) 50.0 (deg) 


high humidity 


direct shear ses ee 51.0 (deg) 
triaxial saturated 2.2 ~ 2.4 (MPa) 2.3 (MPa) 
direct shear | high humidity | 2.35 (MPa) 
room 


high humidity 


direct shear 
room 


The moisture conteni of sandstone has a significant influence on the strength of sandstone. As can 
be seen from Figures 3 and 4. the strength of sandstone reduces as the moisture content increases. 
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actions imposed by engineering works. This behavior is usually anributed, in the scarce literature 
available, only to the expansion process due to hydration and tansformation of anhydrite into 
gypsum. There are no detinite proofs in favor of this theory and it is questioned by some authors. 
The study of rheological properties, in both viscoelastic and viscoplastic fields, may bring some 
light to explain those deformacions. Preliminary results of an experimental program including long- 
term unconfined compression tests and axial free swelling tests on evpsum and anhydrite samples 
from Brariland Argentina are aresanted The results indicate that 4:2 laste behavior of 
these rocks is significant to be considered with respect to the siress increase imposed by some 
ign due to by dration must be considered as 


well. 


1 INTRODUCTION 


Evaporitic rocks present some zeomechanical peculiar szhavior which has called the attention of 
researchers from several areas. Abundant studies on seit rock are an example. In spite of broad 
knowledge about salt rock, there is only rare and sparse inform2cion in the literature abou: the 
geomechanical behavior of sutface evaporitic rocks fixe anhyérze and gypsum, The literature 
contains information about some ancient underground works loczzed in the Alps, which presente! 
long term deformations, large enouzh to cause significant damaze to the reinforced concrete 
(Steiner. 1993; Sahores, 1963: Henke. 1976; Grob, 1975; Krause. 1376). 

Several authors explain this phenomenon based oa the theors of anhydrite hydration (CaSO.) 
transforming it into gypsum (C2SO..2H-0). 

This mineralogical process may involve volume increase of the order of 60%. This theory 
Colides with two important facts. First, a set of important: underground works excavated in this Type 
of rock did not present swelling under that level of deformation (Sahores. 1962). Second. 
geologival evidences of the minevalogical transtormation do not show the occurrence of such 
volume inerease (Orti Cabo an Ortiz, 1981: Carazzi, 1993 

Facing a theory without definite prools, the s 
offer another point of view. T sported here, encompassing 
all the stages of geomechanies! characterization of these rocks. ranging from conventional shert 
duration to long term tests (eres under uniaxial compression ane ¢ 
items, relevant preliminary restits ese tests will be presented. 


ropertizs of these rocks will 


ed 


2 MATERIALS AND METHODS 
Samples of anhydrite and gypsum from the Auquilco Formation of lower Jurassic age from 
Mendoza, Argentina (Auquilco anhydrite, alabastro gypsum, laminar gyosum and massive _ 
gypsum) were colected from the Principal Andean Cordillera. Gypsum from the Santana Formation 
of lower Cretaceous age, Araripe Basin, from Pernambuco, Brazil (nodular and banded gypsum) 
and anhydrite from the Muribeca of Aptian age from Sergipe, Brasil where also colected. 
Petrographic analyses, determination of index properties and short duration tests (Brazilian, 
uniaxial compression and triaxial compression tests) were carried out according to ISRM 
Suggested Methods. The latter were performed in a servo-controlled testing machine MTS 815 
from the Rock Mechanics Laboratory, Sdo Carlos Engineering School, University of Sao Paulo. 
For the free expansion tests, the method suggested by ISRM was adopted. The creep tests were 
performed according to the experimental standards proposed by ASTM (D4405-93). 


3 GEOTECHNICAL PROPERTIES 


Tablel. Main geotechnical characteristics of the rocks studied. 
Table 1. Geotechnical properties 


Tensile Strength Uniaxal Strength Defor. Modulus Poisson Ratio 
Rock MPa MPa GPa - 
EE SE bi Se a ee ee 
Alabastro Gypsum 0,5+0,1 4821 3,3 + 1,27 0,340.26 


25.19 +6,1 19,67 +0,9 
32,5 + 12,3 
30,67 +4,5 
36,7 £52 


69,2 +23 


0,21 +0,01 
0.26 + 0,18 


Banded Gypsum 
Nedular Gypsum 
Auquilco Anhydrite . 


0,2 + 0,04 
0,21 £ 0,02 


6.8225 


102 + 11,2 ° 0.27 + 0,03 


4 TIME DEPENDENT BEHAVIOR 


a1 Free Axial Swelling Tests 


In this research. five axial swelling test were carried out on Auquilco anhydrite samples and four 
on gypsum. The tests consisted of measuring axial displacement of specimens submerged in a 
water solution of catcitim sulphate. Figure 1 shows the results of those tests. The first tests on 
anhydrite (curves | and 2) and gypsum (curves 7, $8 and 9) were carried out by initially distilled 
water, and then satured solution, in order to investigate the effect of concentration on the swelling 
process, 

it can be clearly observed that anhydrite samples (curves | to 5) and gypsum samples with high 
anhydrite content (curve [0) show linear increase of the axial strain and gypsum samples with low 
anhydrite content show volume decrease, probably due to gypsum solution, Anhydrite swelling 
shows a linear trend with time and strain rates vary between 0.6 and 1.1% per year. Macroscopic 
observations of specimens after tests show whitish surface and the formation of gypsum acicular 
crystals forming a continuous layer, Microscopic observations are not available as yet with respect 
to the reaction interface. [t is possible that anhydrite specimen swelling is due in large amount to 
the volume increase from this extemal gypsum layer, and not to generalizated transformation of the 
anhydrite matrix. The transformation seems to occur Only in a thin external layer responsible for 
the whitish color at the surface. On the other hand, gypsum specimens did not show this external 
gypsum layer, and the volume decrease may be'a consequence of gypsum and carbonate solution. 
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4.2 Creep test under constant stress uniaxial compression 


Axial free swelling strain vs time —— Annydrite 


~ +++ + +Gypsum 
Gyp+Anh 


Axial Swelling Strain 
(%) 


0 50 100 150 200 250 300 350 400 450 500 550 600 
Time (days) 


Figure 1. Axial free swelling strain vs time for anhydrite and gypsum 


Creep tests under constant stress uniaxial compression are underway in order to characterize the 
stress-strain-time behavior in a column-type frame which allows five test to be carried out 
‘simultaneously. The frame, the loading system and the data acquisition system are installed in a 
controlled temperature and humidity room (18°C £1°C and 37% £ 5% relative humidity). The test 
consist of three loading and unloading cycles with duration of 22 and 7 days respectively, The 
applied loads correspond to stress ranging from 25% to 80% of the unconfined compressive 
Strength. Viscoelastic and viscoplastic parameters were evaluated for the primary and seconda 
creep condititons. Burger’s rheological model was adopted and the constants were obtained by 
curve adjustment to the data. The curves of axial and radial strain versus time for banded and 
nodular gypsum and Auquilco and Muribeca anhydrites are shown in Figure 2. 

According to Burger's model, axial strain rate is given by: 


-Gr 
PP ae ed Se en 
OK 3G, 3G, 3G; 33 


where K is the time-independent bulk modulus, G; and G; the shear moduli of delayed anc 
instantaneous responses respectively, 7; and 4; are viscosity components controlling primary ané 
secondary creep. Curve adjustement by the least square technique yielded the values shown in 
Table 2 for different stress levels, 

For Muribeca anhydrite, G, varies from 57 GPa for stress level corresponding to 0.4c. to 33+ 
GPa for 0.8q.. indicating that the component of delayed elastic strain decreases with the intens 
of applied stress: 7; varies between 130 and 230 GPa.day forthe same stress level; Gy is roughh 
constant, varying in the range of 23 to 29 GPa ; n: increases linearly with the stress level, Straix 
rates ranged between 6.05x10%days" (7x10"°s") and 8.64x10“days™ (1x10°'s"). For stress lower 
than 20MPa viscous flow is not apparently observed. Plastic strains, as observed in Figure 2. 
increase with time, and viscoelastic strains keep constant, about 10% of total strain. 

Auquilco anhydrite presents time-dependent strain components less significant than Muribeca 
anhydrite, The former has higher strength and lower deformability. Strain rates vary in the range of 
8.6x10%days” (1x 1075") and $.6xl07days" (1x10°s") and total viscoelastic strain in the range of 

“) to 11%, with an average of 3%, Plastic strains are not much significant and are in the range of 
7% (Figure 2). Other important aspect is the fact that it does not present viscous behaviour for 
Siresses higher than 25MPa (0.30.). The difference in geomechanical properties presented by twe 
types of anhydrite is due to ditferent geological and deformation histories. Auguileo anhydrite was 
Subjected more severe tectonic processes stiffening the material. Gypsum presents strain-time 
Curves with well detined primary creep and secondary creep presents low rate, Strain rates for both 
types of gypsum ranges between 1.7x107days” (2x10°s") and 


© 
' 
‘ 
® 
1 
' 


4 
RRRNE 
{Hh 
WE: 
g 
5 
$ 
Zz 
< 
Pa 
Ww 
2 
[eq 
= 
= 
— 
RaaS 
2322 
RExR 
+ ~ { 
is 
Hi | 
GuARS 
322 
i 
t< 
| 


w 
| 
= 
ra] 
= 
2 
< 
fe) 
oO 
= 
5 
S 
3 
z 


396 


Teta 


£ 


vob wy 


a 


4. 


but Oa 
sot 8 
gue ea 

sveae f, 
coco B 


4oL ut 


netod t 


TIME (DAYS) 


Tite (DAYS) 


| 


= PARAS (03.2) —e— GIPARBAS (18,6) | 


—e— SWAKHAI (12.2) —92—GIPARDAS (38) 


GIPARMAI (13) —8— GIPARMAI (14.6) 


| —ne— GIPAFMA 1 (21,6) —@— GRARMAZ (13) 


Ts 


NODULAR GYPSUM 


b= GIPAMBA?2 (8.4) 
O GIPAREAN (9.3) 


GIPARHA? (12,6) 


oe GIPAIAT (N24) 
x 


BANDED GYPSUM 


[4 GPAAZ (215) —w—GIPARMA2 (25.7) 


0.00120 


ooo1o0 


0.00020 
0.00000 


: 


#oous0 
eg0u00 
asuoag 
uguuze F 


-0,00020 


aude 


0.00040 


ome 


enovey | 


TIME (DAYS) 


TIME (DAYS) 


- time curves of the anhydrite and gypsum 


Figure 2. Deformation 
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1.7x10%days' (2x10%s"). The ternary diagram of Figure 3, modified from Ozgenoglu (1998), 
shows that the viscoelastic strain component ranges between 5% and 20% for banded gypsum and 
between 8% and 17% for nodular gypsum. 


Table 2. Burger's Model Parameter 
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GYPSUM 


Figure 3. Ternary diagrams for anhydrite and gypsum 


Plastic strain components are in the range of 17% to 23s for the banded variety and $%o to 2574 for 
the nodular one. Both types of rock do not seem to present a threshold stress level for the onset of 
creep. 

[It was observed by direct and indirect means that the initiation crack propagation is the 
mechanism responsible for creep strain. [t is plausible to state that a constitutive model based on 
the theory of progressive damage will be adequate to predict stress-strain-time behaviour of t 
type of rack, 


3 FINAL REMARKS 


1. The process of transformation of anhydrite into gypse not clear. and depends. in real 


present volume increase by surface deposition of gypsum crystal by mineralogical conversion. 
Experimental results show that the formation of this gypsum layer only occurs for anhydrite and 
not for gypsum. Strain rates range between 0.6% to 1.1% per year. For a !0-m diameter tunnel and 
swelling zone of 2 m around the tunnel, convergence would be in the range of | to 2 cm per year. 
These values are in the range of what has been reported by Steiner (1993) for the Hauenstein 
Summit tunnel (>10mm‘year) and Wagenburg Tunnel (7mm/year). 

2- Creep tests indicate that gypsum, as expected, present higher viscoelastic and viscoplastic 
strains than anhydrites. Strain rates during secondary creep are roughly 10 times larger for gypsum 
as compared to anhydrite. Strain rates for Muribeca anhydrite (10° to 10°'°s"') are larger than these 
for Auquilco anhydrite (107 to 10s"). This difference can be attributed to the higher level of 
stress history to which the latter has been subjected. Considering free swelling rates for Auquilco 
anhydrite and the same conditions for the hypothetical tunnel, convergence rate in the range of 0.6 
to 6 cm per year could be expected. 

3- Comments of item 2 and 3 show that free swelling and creep strain are of the same order of 
magnitude. The large strains observed in tunnels can be attributed to them, ocurring isolately or in 
combination. 

4- Rheological models are a simple form to obtain deformation parameters. However they fail 
for generalized conditions like 3-dimensional stress field. Some of them, like Burger's model, fails 
even for nonlinear behaviour under uniaxial stress. Future analyses should focus on the 
microcracks propagation phenomenon for the proposition of a constitutive model based on physical 
mechanism. 

5- Since both creep and transformation of anhydrite into gypsum are phenomena that can occur 
simultaneously in nature, it seem to be appropriated to study experimentally coupled effects 
between the nwo processes. 
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Statistical evaluation of rock mass properties at TBM driving 
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ABSTRACT: The accurate and adequate input data obtained by engineering-geological 
investigation are necessary for appropriate evaluation of the rock environment. Results of common 
methods of engineering-geological investigation in real tunnelling conditions do not provide 
sufficient number of the measurement results for a proper use of statistical methods, and the results 
are only available with a certain time delay. This disadvantage may be eliminated by the inverse 
continuous method for rock mass properties determination (IKONA). The model rock mass 
strength obtained by this method continuously simulates the compressive strength changes in the 
tunnel face. The paper analyses the model strength by statistical methods for the purpose of rock 
mechanies knéwledge specification The paper presents the results of the statistical evaluation of the 
Branisko rock mass properties, where the first motorway tunnel in Slovakia has been driven by the 
tunnel boring machine (TBM). 


| INTRODUCTION 


The base of rock mechanics analysis is the geological database, which covers rock types, structural 
discontinuities and rock material properties. Proper rock environment assessment results require the 
precise and adequate input values obtained usually by engineering-geological investigation. 

During full profile machine tunnelling there is a possibility to record and evaluate the tunne! 
driving process parameters affected by rock mass properties. This full-profile tunnelling process 
monitoring provides the observation of the rock properties changes on the tunnel face, thus 
complementing the engineering-geological investigation. 

During the excavation of easter part of the Branisko exploratory gallery, a unique computer 
monitoring and optimising system. called WORS, was installed on the Wirth TBII-330-H™M 
tunnelling machine. The principle of the WORS system has been developed at the Institute of 
Geotechnies, Slovak Academy of Sciences, and has been deseribed in several works of Krupa 
(1998). In brief outline, the WORS system enables optimisation of the TBM operation by scanning 
the TBM parameters (applied thrust. torque, revolutions and position of cutterhead) with following 
real-time calculations of net advance rate, specitic disintegration energy and recommended thrust. 
Driving with recommended thrust guarantees the excavation with minimum specific volume 
disintegration energy. 

Another nonpareil project, developed at the Institute of Geotechnies, SAS, is the [KONA 
system. Using the investigation of the rock-tool interaction, a mathematical model of the 
interaction has been designed by Krupa (1998), giving opportunity for determination of rock mass 
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